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Summary 


Consider  a  one-on-one  engagement  between  two  adversarial  combat  platforms,  and  an  important 
outcome  of  that  engagement:  who  is  victorious  and  who  is  not.  A  simple  answer  is  that  the 
victor  lives  while  the  vanquished  dies,  but  for  a  meaningful  analysis  of  platform  survivability, 
this  answer  is  insufficient.  A  better  answer  ought  to  reflect  a  more  practical  view  of  victory  and 
defeat;  one  that  illuminates  the  fate  of  each  platform  if  we  could  repeatedly  conduct  the  same 
engagement.  Given  that  the  current  paradigm  for  platform  survivability  is  one  of  layered 
survivability,  there  is  a  need  for  a  unified  framework  within  which  one  can  analyze  the  full  range 
of  survivability  options  layered  survivability  suggests.  This  paper  sketches  such  a  framework, 
one  built  upon  a  mathematical  construct  called  the  stochastic  duel.  As  we  construct  this  duel,  we 
first  consider  a  platform  that  engages  a  passive  target,  and  we  model  this  process  via  the 
techniques  of  renewal  theory.  From  there,  we  model  the  one-on-one  engagement  by  equating  it 
to  a  series  system,  a  choice  that  allows  us  to  apply  the  method  of  competing  risks.  Within  the 
competing  risk  framework,  we  then  derive  a  metric  for  a  platform’s  probability  of  survival. 
Finally,  we  extend  this  formalism  to  include  survivability  measures  by  relying  on  the  concept  of 
layered  survivability. 
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1.  I  ntroduction 


This  paper  considers  a  simple  one-on-one  engagement  between  two  combat  platforms  in  an 
effort  to  study  the  outcome:  who  is  victorious  and  who  is  not.  We  find  the  trivial  answer,  that 
the  victor  lives  while  the  vanquished  dies,  insufficient  for  our  analytical  needs.  The  answer  we 
seek  ought  to  reflect  a  more  practical  view  of  victory  and  defeat;  one  that  illuminates  the  fate  of 
each  platform  if  we  repeatedly  conducted  the  same  engagement.  We  shall  assume  that  the  victor 
moves  on  to  other  engagements  and  the  vanquished  does  not;  however,  we  will  not  assume  the 
destruction  of  the  vanquished  or  that  the  victor  proceeds  unscathed.  In  other  words,  we  assume 
that  the  victor  remains  mission  capable  while  the  vanquished  does  not.  The  U.S.  Department  of 
Defense  Dictionary  (2001)  defines  mission  capable  as  the  “material  condition  of  an  aircraft 
indicating  that  it  can  perform  at  least  one  and  potentially  all  of  its  designated  missions.”  While 
this  definition  applies  to  aircraft,  we  take  this  definition  to  apply  equally  to  ground  combat 
platforms.  Accordingly,  how  can  one  determine  each  platform’s  fate  if  we  could  conduct  this 
single  engagement  repeatedly?  We  motivate  this  research  with  the  desire  to  address  this  very 
question;  thus,  we  turn  to  a  mathematical  model  called  the  stochastic  duel. 

Stochastic  duels,  as  an  analytical  approach,  were  first  proposed  by  Williams  and  Ancker  (1963). 
Their  approach  involved  two  duelists  who  fired  at  random  intervals.  Each  duelist  possessed 
unlimited  ammunition,  fired  with  a  fixed,  single  shot  probability  of  hit,  and  the  first  one  to  hit  the 
other  became  the  victor.  Williams  and  Ancker  used  the  term  marksman  problem  to  describe  a 
combatant  firing  at  a  passive  target  until  he  observes  a  hit.  Then,  by  assuming  stochastic 
independence  between  these  marksmen,  Williams  and  Ancker  turned  these  marksmen  against 
one  another  to  form  what  they  termed  the  fundamental  duel. 

We  set  two  goals  for  this  paper:  (1)  explore  the  fundamental  duel  and  the  conditions  that  impose 
stochastic  independence,  and  (2)  incorporate  survivability  measures  into  the  duel.  At  the 
conclusion  of  this  paper,  we  will  have  built  the  foundation  for  a  follow  on  paper  that  describes 
our  model  of  the  fundamental  duel.  We  structure  this  paper  first  to  express  stochastic  duels  as 
renewal  processes,  next  to  study  the  fundamental  duel,  and  finally  to  add  survivability  measures 
to  the  renewal  duel  formulation.  Consequently,  section  2  provides  a  brief  review  of  renewal 
theory,  and  in  section  3,  we  employ  renewal  theory  to  model  the  component  processes  of 
stochastic  duels  as  renewal  processes.  Section  4  constructs  the  fundamental  stochastic  duel  of 
Williams  and  Ancker  by  framing  the  duel  as  a  competing  risk.  In  section  5,  we  employ  the 
mathematics  developed  in  the  previous  sections  to  derive  solutions  for  the  fundamental  duel; 
solutions  that  we  then  validate  by  comparing  them  to  results  obtained  by  Williams  and  Ancker. 
Next,  section  6  introduces  a  concept  called  layered  survivability  in  order  to  establish  a  basis  for 
section  7,  which  extends  our  renewal  duel  formulation  to  include  survivability  measures. 
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2.  A  Brief  View  of  Renewal  Theory 


We  begin  our  review  by  considering  a  sequence  of  arbitrary  events  that  occur  at  random 
intervals.  We  assume  that  one  can  model  these  intervals  as  independent  random  variables  with 
identical  distributions.  When  the  length  of  the  interval  between  successive  events  is  non¬ 
negative  with  probability  one  (i.e.,  almost  surely),  we  term  the  sequence  a  renewal  process  (Kao, 
1997).  As  we  describe  renewal  processes,  we  rely  on  the  standard  example  of  that  process:  a 
sequence  of  light  bulbs  where  each  bulb  runs  consecutively.  At  time  S0  =  0,  we  apply  power  to 
the  first  bulb.  After  a  period  of  operation,  the  light  bulb  fails  and  we  let  St  record  the  epoch  of 
first  failure.  We  immediately  replace  the  bulb  with  an  equivalent,  and  replicate  this  process 
when  the  bulb  fails  at  S2,  S3,  etc.  Our  desire  is  two-fold;  first,  we  seek  a  count  of  replacements 
up  to  some  time  t,  and  second,  we  wish  to  understand  how  the  process  behaves  as  t  grows 
without  bound.  We  set  S0=  0  and  let  [  Si  }”=1  denote  the  times,  or  equivalently  epochs,  at  which 
the  events  occur.  We  define  the  interval  between  successive  events  to  be  X(  =  St  -  S(_t ,  and 
assume  that  the  lengths  of  these  intervals  are  independent  and  identically  distributed  (i.i.d.) 
random  variables  with  F(x )  =  P(Xt  <  x).  Thus 

a) 

1=1 

is  the  elapsed  time  until  the  occurrence  of  the  nth  event.  Set 

N (t)  =  max{n :  S n  <t]  (2) 

so  that  N(t)  counts  the  events  in  the  interval  (0,  t ].  Thus,  we  refer  to  {N(t),t  >  0}  as  the 
counting  process,  and  {Sn,n  >  0}  as  the  partial  sum  process',  however,  in  keeping  with  common 
practice,  either  process  may  be  referred  to  as  the  renewal  process  (Karlin  and  Taylor,  1975;  Kao, 
1997).  We  may  determine  the  distribution  of  N(t)  from  our  knowledge  of  F(x)  by  employing 
the  following  relationship  given  by  Karlin  and  Taylor 

N(t)>n  iff  Sn<t,  (3) 

where  iff means  if  and  only  if.  An  immediate  consequence  of  this  relationship  is  that 

P(N(t)>n)  =  P(Sn<t) 

=  Fn(t),  t  >0,n  =  l,2,---, 

where 
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[ F  At- r)dF (r),  n>  1, 

Fn(t)  =  \{  -lV  (5) 

Fit),  n=  1. 

Consequently,  Kao  (1997)  determines  that 

'P(N(t)>n)-P(N(t)>n  + 1), 

P((V<0  =  «)=P(S,<()-i,(S, „,<().  (« 

P,(0-P„1«>.  (>0,n  =  l,2,'". 

We  now  turn  to  the  renewal  function  M(f),  which  we  define  as 

M(t)  =  E[N(t)].  (7) 

By  standard  arguments,  various  authors  rewrite  eq  7  as 

oo 

M(t)  =  YjnP(N(t)  =  n) 

n= 1 

00 

=  2>(iV(f)>n)  (8) 


00 

=2X«>. 

/;=1 

Karlin  and  Taylor  (1975)  formally  prove  convergence  of  the  series 

oo 

Mit)  =  Y,Fnit).  (9) 

n= 1 


We  can  examine  the  renewal  function  in  more  detail  by  recognizing  that  the  underlying  counting 
process  is  a  probabilistic  replica  of  itself  between  each  event.  If  we  denote  the  first  epoch  as  Xv 
and  condition  on  its  occurrence  at  X,  =  x,  we  can  write 


£[W(OIX1=a:] 


jl  +  M(t-x),  x<t, 
[O,  x  >  t, 


(10) 


where  the  quantity  1  +  M(t  -x)  counts  the  first  event  plus  any  remaining  events  in  the  interval 
(at,  t\.  Thus,  by  applying  the  law  of  total  probability  (Bartoszynski  and  Niewiadomska-Bugaj, 
1996)  over  the  range  of  arrival  times,  we  find  a  functional  equivalent  for  the  renewal  function 
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M(t)  =  E[N(t)] 

oo 

=  J^[iV  (01  =  x]dF(x) 

0 

=  J[1  +  M(t-x)]dF(x)  (U) 

0 

t 

=  F(t)  +  jM(t-x)dF(x). 

o 

We  will  find  a  solution  to  eq  1 1  after  considering  Theorem  1. 

Theorem  1 

Suppose  a  is  a  bounded  function.  There  exists  one  and  only  one  function  A  bounded  on 
finite  intervals  that  satisfies 

t 

A(t)  =  a(t)  +  $A(t-T)dF(T).  (12) 

0 


This  function  is 


A(t)  =  a(t)  +  \a(t-T)dM(T).  (13) 

0 

where  M(t)  =  Fn  (t)  is  the  renewal  function.  (Karlin  and  Taylor,  1975,  p.  184-185). 

Observe  that  when  a(t)  =  F(t )  and  A(t)  =  M(t),  eq  1 1  satisfies  the  requirements  of  Theorem  1 . 
Thus,  we  can  write  the  following  unique  (and  equivalent  by  the  property  of  convolutions) 
solutions  to  eq  1 1 

t 

M(t)  =  F(t)+  \F(t-x)dM(x) 

:  (i4) 

=  F(t)  +  ^M(t-x)dF(x). 

o 

Equation  14  completes  our  review  of  renewal  theory,  and  we  finish  the  section  by  noting  that  in 
general,  closed-form  solutions  to  this  equation  are  only  possible  with  the  simplest  of  inter-arrival 
distributions.  We  now  begin  our  study  of  stochastic  duels  by  employing  renewal  theory  to 
mathematically  express  the  component  processes  of  the  fundamental  duel. 
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3.  Stochastic  Duels:  The C omponents as  Renewal  Processes 


In  this  section,  we  study  a  one-on-one  engagement  between  two  adversaries;  one  arbitrarily 
termed  the  combatant  and  the  other  the  target  simply  to  signify  that  the  target  is  the  object  of  the 
combatant’s  attention.  The  combatant  defeats  the  target  by  undertaking  and  successfully 
completing  three  consecutive,  but  distinctly  different,  processes.  To  describe  these  processes,  we 
will  introduce  some  terms  as  convenient  referents  for  the  mathematics.  In  the  first  of  these 
processes,  the  combatant  seeks  out  a  target;  we  incorporate  this  process  into  the  generalized 
marksman  problem.  Here,  we  take  up  the  term  “marksman  problem”  in  the  sense  of  Williams 
and  Ancker  (1963),  but  adapt  it  to  reflect  our  extension  to  their  work.  Secondly,  once  the 
combatant  detects  a  target,  he  then  engages  by  loading,  aiming,  and  firing  his  weapon  at  the 
target;  we  call  this  process  the  marksman  problem  to  reflect  that  this  process  represents  the 
problem  as  originally  conceived  by  Williams  and  Ancker.  Finally,  once  the  combatant  hits  the 
target  we  want  to  ascertain  the  extent  to  which  the  target  sustains  damage.  We  employ  the  term 
damage  process  to  describe  the  accretion  of  damage  by  a  target  with  each  hit. 

In  this  section,  we  will  first  express  the  marksman  problem  as  a  renewal  process.  Then,  we  will 
generalize  this  model  to  include  the  process  of  detection  by  considering  the  marksman  problem 
as  a  delayed  renewal  process.  Finally,  we  will  model  the  damage  process  via  renewal  theory. 
Thus,  after  finishing  this  section,  we  will  obtain  a  complete  mathematical  model  for  one-half  of  a 
one-on-one  duel  by  joining  a  generalized  marksman  problem  and  a  damage  process  together. 

3.1  T  he  M  ar ksman  Problem 

We  begin  by  considering  a  combatant  as  he  fires  at  the  strictly  passive  target.  The  combatant 
begins  a  trial  at  time  zero  with  weapons  unloaded.  After  a  random  interval,  the  combatant  loads 
and  fires  at  the  target.  Once  the  combatant  fires,  he  instantaneously  observes  the  outcome  of  his 
shot.  If  he  misses  the  target,  he  repeats  the  firing  process;  however,  if  he  hits  the  target,  he 
begins  a  new  trial.  Figure  1  provides  a  sample  realization  of  this  process  where  we  employ 
dashed  arrows  to  signify  shots  that  miss,  and  solid  arrows  to  signify  shots  that  hit.  Figure  2 
represents  the  associated  counting  process  for  the  hits  within  this  sample  process.  If  our  sole 
interest  lay  in  the  distribution  of  the  intervals  between  consecutive  shots  taken  in  (0,f],  we  could 
avail  ourselves  of  the  approach  outlined  in  section  2  to  determine  that  distribution.  However,  in 
our  case,  the  outcome  of  each  shot  is  also  a  random  process.  Our  interest,  therefore,  ought  to  be 
in  deducing  the  distribution  of  intervals  between  consecutive  occurrences  of  the  same  outcome. 
Before  we  can  find  that  distribution,  we  first  make  some  observations  regarding  this  process  and 
then  define  some  necessary  variables. 


6 


Figure  1.  Sample  shot  history. 


One  could  view  the  combatant’s  sequence  of  shots  and  shot  outcomes  as  a  stochastic  process  of 
the  form  {S‘,Y.},  where  5  s  denotes  the  firing  time  of  the  shot,  Y  denotes  the  location  of  the 
shot  impact  (not  necessarily  on  the  target),  and  i  indexes  the  shot.  While  perhaps  desirable, 
modeling  such  a  process  is  practically  impossible  for  all  but  the  most  simplistic  functions. 
However,  if  we  model  this  process  as  { S- ,  /  (F(. ) } ,  we  can  resort  to  established  techniques  to  find 
our  solutions.  Here  we  employ  the  indicator  function  /(¥,.)  in  the  following  sense 

/(F)  {*’  ^  the  point  of  shot  impact  is  on  the  target, 

[0,  if  the  point  of  shot  impact  is  off  the  target. 

By  employing  the  indicator  notation,  we  recognize  that  we  are  dealing  with  a  stochastic  process 
that  consists  of:  (1)  a  shot  process  with  an  embedded  (2)  hit  process.  In  this  process,  shot  events 
occur  at  epochs  {5/  }Jl,  and  hit  events  occur  at  epochs  {S*  }*=1.  We  note  that  for  every  5*,  there 
is  one,  and  only  one  S- ,  such  that  S*  =  5 / ;  however,  in  general  the  converse  is  not  true. 
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We  start  our  definitions  by  setting  S*  =  5*  =  0  so  that  we  can  denote  the  interval  between 
consecutive  shots  as  X.  =  S ■  -  Sj_ and  define  X'j.  =  S%  -  to  be  the  interval  between 
consecutive  hits.  Figure  3  depicts  the  relationship  between  these  variables  as  they  apply  to  the 
sample  process  of  figure  1.  Consequently,  if  each  X-  >  0  almost  surely,  we  can  model  this 
process  as  a  renewal  process.  We  will  assume  that  the  intervals  X.  and  outcomes  F  form  an 
i.i.d.  sequence  of  random  variables  of  finite  variance,  with  distributions  Fx  (x)  =  P(X-  <  x)  and 
Fy  (y)  =  P(Yi  <  y)  respectively.  We  note  that  while  Karlin  and  Taylor  (1975)  allow  for  the 
possibility  that  F  may  depend  on  X. ,  we  shall  assume  independence  between  F  and  X. ,  and 
defer  our  exploration  of  dependence. 


Figure  3.  Relationships  among  key  variables. 

By  employing  the  indicator  function  as  we  did,  we  restricted  the  possible  outcomes  for  a  given 
shot  to  two:  either  the  shot  hit  or  it  missed.  Thus,  we  can  employ  Doob  (1994)  to  define  the 
probability  of  hit  as 

ph=  f  dFY(y).  (15) 

/(h)=i 

While  literature,  such  as  the  work  of  Helgert  (1971),  Jaiswal  (1997)  or  Przemieniecki  (2000), 
exists  on  computing  the  probability  of  hit  for  many  types  of  targets,  for  the  purpose  of  this  paper, 
we  shall  assume  that  the  probability  of  hit  is  both  known  and  constant. 

With  our  definitions  complete,  we  can  now  develop  the  marksman  problem  as  a  renewal  process. 
Our  goal  is  to  find  the  probability  distribution  that  describes  the  interval  between  consecutive 
hits,  a  goal  we  represent  by  constructing  the  renewal  argument  M(t )  =  P{X\  <  t ).  Our  task  is  to 
find  a  functional  equivalent  for  this  argument  based  upon  our  knowledge  of  Fsx  (x).  We  begin 
by  conditioning  on  the  epoch  of  the  first  shot 
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(16) 


[M{t)  \X‘ 


|{P(X1,‘<OI/(F1),X1s=x}, 

1°. 


x  <t, 
x  >t. 


In  eq  16,  the  term  { P(X *  <  t)  I  J(Fj),  X[  =  x}  reflects  that  our  argument  necessarily  depends  on 
the  outcome  of  the  first  shot.  As  noted  before,  we  restricted  the  range  of  possible  outcomes  for  a 
given  shot  to  two:  either  a  hit,  which  implies  the  event  I(Yl)  =  1;  or  a  miss,  which  implies  the 
event  I(YX)  =  0.  Clearly,  if  the  first  shot  hits  the  target,  then  P(X f  <t)  =  \  and  we  restart  the 
process.  However,  if  the  first  shot  misses  the  target,  then  P(X'‘  <  t )  becomes  the  probability 
that  a  hit  occurs  in  (x,f] ,  or  simply  P(X[l  <t)~  P(Xf  <  x),  which  is  equivalent  to  M(t  -  x). 
Thus,  after  conditioning  on  the  outcome  of  the  first  shot,  we  write 


{P(X?  <01/(7,),^ 


|  M(t-x), 


I(Yl)  =  0, 

/(^)  =  1. 


(17) 


We  then  apply  the  law  of  total  probability  over  the  distribution  of  7(Fj) ,  and  combine  this  result 
with  eq  15  to  find 


[M(0  I  X[  =  jc]  =  { {P(Xf  <  t )  I  I(Yx)}dFI(Y) 
=  Ph  +  0--Ph  )M(t-x). 


(18) 


By  combining  eqs  16  and  18,  we  can  now  describe  the  state  of  our  process  after  the  first  shot 


[M(t)  \X[ 


]  =  |  Ph+Q-PhWV-x), 


x  <t, 
x  >t. 


(19) 


Again,  by  applying  the  law  of  total  probability,  but  this  time  over  a  range  of  inter-shot  arrival 
times,  we  find  our  functional  equivalent  for  the  renewal  argument  M(t)  =  P(X*  <  t)  in 

oo 

M(t)  =  J  [_M(t )  I  X[  =  x]dFsx  (jc) 

0 

t 

=  J  [ph  +  (l-ph)M(t-x)]dF’(x)  (20) 

0 

t 

=  phFsx  (t)  +  {l-ph)\M(t-  x)dFsx  (x). 

o 


We  leave  eq  20  as  our  mathematical  description  of  the  marksman  problem,  and  the  principal 
result  of  this  section. 
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3.2  G eneralizing  the  M  arksman  Problem 

We  proceed  as  in  section  3.1  with  a  combatant  and  the  combatant’s  unloaded  weapon  starting  a 
trial  at  time  zero;  however,  now  the  combatant  must  first  search  the  battlefield  to  find  the  target. 
We  shall  employ  the  same  notation  as  the  previous  section  by  letting  {5/,/(Ff)}  be  our 
stochastic  process.  Like  the  marksman  problem,  we  will  model  this  process  as  a  renewal  process 
in  which  shots  occur  at  epochs  {5/  }"p  and  hits  occur  at  epochs  {5*  }*=1.  We  will  let 
X. ,  X%,  and  F  take  the  same  meanings  as  in  section  3.1.  F  will  have  the  same  distribution 
FY  (y)  =  P(Yj  <  y ),  as  will  X.  for  i  >  1  in  the  distribution  Fx  (at)  =  P(X.  <  x).  However, 
unlike  the  marksman  problem,  the  first  epoch  now  has  the  distribution 

Gsx(x)  =  P(Ssl<x)  =  P(  F'  <  x),  and  we  assume  this  time  represents  the  time  to  first  detect  and 
engage  the  target.  Karlin  and  Taylor  (1975),  Kao  (1997),  and  others  term  the  class  of  processes 
represented  by  this  phenomenon  delayed  renewal  processes. 

Kao  (1997)  gives  the  renewal  function  for  a  delayed  renewal  process  as 

oo 

MD(t)  =  E[ND(t )]  =  £G(0*iv,(0.  (21) 

n= 1 

In  eq  21,  we  employ  the  subscript!)  to  signify  delayed  renewal  formulation;  ND(t )  denotes  the 
delayed  counting  process;  denotes  convolution;  G  and  F  are,  respectively,  the  distributions  of 
the  first  and  subsequent  inter-arrival  times;  and  the  subscript  n- 1  denotes  the 
(n-l)-fold  convolution  of  F  (see  eq  5).  We  examine  this  process  by  starting  with 
eq  21 

oo 

n-l 

00 

=g(o*2X(o 

M=°  .  (22) 

=g(o+g(o*2X(o 

m=  1 
t 

=  G(t)  +  $M(t-x)dG(x). 

0 


where  M(t)  =  ^Fjt).  Equation  22  shows  that  MD(t )  is  the  solution  of  a  renewal  equation 
of  the  form 

t 

MD(t)  =  G(t)  +  \MD(t-x)dF(x).  (23) 

o 

One  can  further  show  that  both  MD(t )  and  M(t)  obey  Theorem  1  (Karlin  and  Taylor,  1975, 
pp.  198-199),  and,  in  fact,  MD{t)  =  M(t). 
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We  let  the  goal  for  this  section  be  the  desire  to  incorporate  the  process  of  target  search  and 
detection  into  the  marksman  problem;  we  are  now  able  to  reach  that  goal.  First,  recall  that 
Gx  it)  and  Fsx  it)  are  the  respective  distributions  for  the  time  to  fire  the  first  and  subsequent 
shots.  Now,  invoke  the  renewal  argument  M(t)  =  P(X%  <  t)  to  model  the  time  between 
successive  hits.  First,  let  eq  23  take  on  the  quantities  G(t)  =  Gsx  ( t ),  Fit)  =  Fx it)  and 
Md  it)  =  M(t).  Next,  compare  eq  23  to  eq  20.  Thus,  we  immediately  find  a  functional 
equivalent  for  our  renewal  argument  in 

t 

M{t )  =  phGx  it)  +  (1  -ph)\M{t-  x)dFx  (x).  (24) 

0 

We  leave  eq  24  as  the  mathematical  description  of  the  marksman  problem  with  detection  and  the 
principal  result  of  this  section.  However,  we  now  also  recognize  that  the  marksman  problem  of 
eq  20  is  simply  a  delayed  renewal  process  (as  represented  by  eq  24)  with  Git)  =  F{t).  We 
therefore  conclude  that  the  delayed  renewal  process  naturally  incorporates  the  marksman 
problem  of  Williams  and  Ancker;  hence,  we  shall  refer  to  the  renewal  model  in  eq  24  as  the 
generalized  marksman  problem. 

3.3  A  Simple  Damage  Model 

In  this  section,  we  will  again  consider  a  combatant’s  sequence  of  shots  as  well  as  their  outcomes. 
However,  unlike  previous  sections,  we  consider  the  progression  of  damage  inflicted  by  shots  that 
hit  the  target.  As  before,  we  have  a  sequence  of  shots  coupled  with  outcomes,  which,  if  one 
ignores  shots  that  miss,  one  can  view  as  a  sequence  of  hits  with  associated  damage  outcomes.  A 
possible  model  for  such  a  sequence  is  a  stochastic  process  of  the  form  {S*,  Ak },  within  which 
Sh  would  possess  a  distribution  described  by  the  generalized  marksman  problem  (eq  24),  and  A 
would  denote  the  damage  caused  by  the  hit.  While  it  is  possible  to  view  the  sequence  of  shots  in 
such  a  manner,  and  sophisticated  software  codes  (BRL-CAD,  2000;  MuVES,  2000;  ORCA, 

2000)  exist  with  which  to  model  the  damage  a  target  sustains  from  a  given  shot,  we  desire  a 
simpler  analytical  model.  Hence,  we  will  view  our  process  as  { Sk ,  I (Ak ) }  and  turn  to  renewal 
theory  for  our  model. 

To  begin  our  investigation  of  damage  effects,  let  us  assume  that  a  given  hit  either  renders  the 
target  mission  incapable  or  causes  no  damage.  One  might  think  of  this  in  terms  of  armor,  in  that 
armor  either  protects  absolutely  or  has  no  effect.  We  model  this  situation  via  the  indicator 
function  I(Ak)  applied  in  the  following  sense 

f  1,  if  the  hit  caused  lethal  damage, 

/(A,)  = 

[0,  if  the  hit  did  not  cause  lethal  damage. 

Our  reliance  upon  the  indicator  notation  reveals  a  stochastic  process  consisting  of:  (1)  a  damage 
process  embedded  within  a  (2)  hit  process.  Continuing  with  the  notation  of  our  previous 
sections,  we  let  the  hit  events  occur  at  epochs  { Sk  }“=1,  and  define  the  interval  between 
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consecutive  hits  to  be  Xk  =  Sk  -  Sk_v  We  introduce  the  sequence  { Sf  }"=1,  and  define  it  to  be 
the  epochs  of  hits  that  render  the  target  mission  incapable.  We  will  denote  the  interval  between 
consecutive  lethal  hits  by  Xf  =  Sf  -  S?_v  Figure  4  relates  these  variables  to  one  another  by 
employing  the  sample  process  of  figure  1  save  only  that  now  we  have  removed  shots  that  missed 
the  target  and  placed  a  diamond  on  those  hits  causing  lethal  damage.  Once  more,  we  note  that 
for  every  Sf  there  is  one,  and  only  one  Sk  ,  such  that  Sf  =  Sk. 


Figure  4.  Relationships  between  key  damage  model  variables. 

We  assume  the  intervals  Xk  and  damage  outcomes  Ak  to  be  mutually  independent  random 
variables  with  distributions  FA  (t )  =  P(Xk  <  t)  and  Fa  (a)  =  P(Ak  <a),  respectively.  Finally, 
we  turn  to  Doob  (1994)  and  define  the  probability  of  lethal  damage  (or  equivalently  the 
probability  of  kill)  as 

Pk=  J  dFA(a).  (25) 

I(At)= 1 


In  addition,  we  will  assume  pk  to  be  known  and  constant  for  this  paper. 


Our  goal  for  this  section  is  to  determine  the  probability  distribution  of  the  interval  between 
consecutive  lethal  hits,  and  we  model  this  goal  via  the  renewal  argument  M(t)  =  P(X‘l  <  t ).  We 
begin  by  conditioning  on  the  time  of  the  first  hit 


[m(0  iz* 


j{P(Z?<OI/(A1),Z1'‘=x}, 

[0, 


x  <t, 
x  >t. 


(26) 


Again,  {P(Zf  <  t)  I  /(At),  Z*  =  x}  reflects  that  the  process  depends  on  the  damage  caused  by 
the  first  hit.  If  the  first  hit  causes  lethal  damage,  then  clearly  P(Xf  <t)  =  1,  and  the  process 
restarts.  Conversely,  if  the  first  hit  does  not  cause  damage,  then  P(Zf  <  t)  becomes  the 
probability  that  a  lethal  hit  occurs  in  (x,  t],  which  is  simply  P(X‘‘  <t)~  P(Zf  <  x),  or 
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M(t  —  x).  Thus,  after  conditioning  on  the  outcome  of  the  first  hit,  we  write 


{P(x?<t)\i(Al),x';  =  x}  = 


I  (A  i)  =  0, 

/(A)  =  i- 


(27) 


We  now  apply  the  law  of  total  probability,  taken  over  the  distribution  of  /(A,  ),  to  eq  27 
[Mit)  I  X*  =  x]  =  {  {P(Xf  <  t)  I  /(Aj),  Z?  =  x  W(A) 


(28) 


By  combining  eqs  26  and  28,  we  can  now  describe  the  process  after  the  first  shot 


[M(t)\Xhi  =  x]  = 


fft  +(i-p*)M(r-x),  x<f, 

1 0,  x  >  f . 


(29) 


Again,  by  applying  the  law  of  total  probability,  but  now  over  range  of  inter-shot  arrival  times,  we 
find 


M(t)  =  |  [Af(0 1  Z*  =  x]  dFx(x) 

0 

t 

=  {  \pk+(\-pk)M(t-x)]dFhx{x)  (30) 

0 

t 

=  pkFhx  (t)  +  (l-pk)\M(t- x)  dFhx  (x). 

0 

We  leave  eq  30  as  the  primary  result  of  this  section,  and  close  by  noting  that  Fx  it)  may  either 
be  a  given  or  obtained  from  the  generalized  marksman  problem  modeled  by  eq  24.  However,  we 
should  also  note  that  eq  30  represents  a  rudimentary,  though  necessary,  first  step  towards 
incorporating  the  effects  of  damage  upon  a  target.  A  more  detailed,  though  still  elementary, 
damage  model  might  be  a  shock  model  where  system  failure  requires  a  succession  of  “shocks”  to 
a  system  to  occur.  Both  Finkelstein  (1996)  and  Gut  (2001)  discuss  various  random  shock 
models.  Finkelstein  bases  his  model  upon  renewal  theory,  and  Gut  allows  for  mixed  shocks 
where  several  small,  or  one  big  shock,  can  cause  failure.  While  either  shock  model  would  likely 
serve  as  the  next  step  to  a  more  detailed  damage  model,  for  the  remainder  of  this  paper,  we  shall 
only  consider  our  rudimentary  damage  model. 
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4.  T  he  F  undamental  D  uel  as  a  C  ompeting  R  isk  Problem 


In  the  previous  section,  we  focused  on  a  single  combatant  as  he  engaged  a  strictly  passive  target 
in  order  to  model  the  components  of  a  stochastic  duel  with  renewal  theory.  Now,  we  will  take 
advantage  of  this  work  to  model  th e,  fundamental  stochastic  duel  of  Williams  and  Ancker  (1963) 
as  a  one-on-one  combat  between  two  marksmen  we  call  Red  and  Blue,  a  process  depicted  in 
figure  5.  We  form  the  fundamental  duel  by  assuming  Red  and  Blue  continue  combat  until  one  is 
no  longer  mission  capable.  This  underlying  assumption  presents  us  a  hidden  benefit  in  our  study 
of  that  duel:  the  ability  to  draw  a  parallel  between  the  fundamental  duel  and  a  “series  system.”  If 
we  equate  a  combatant  that  is  no  longer  mission  capable  to  a  failed  component,  we  realize  that  a 
series  system  that  fails  when  a  component  fails  models  a  combat  that  terminates  when  a 
combatant  becomes  mission  incapable.  When  we  view  a  system,  and  therefore  the  fundamental 
duel,  in  this  manner,  we  recognize  a  specific  instance  of  a  general  class  of  problems  that  David 
and  Moeschberger  (1978)  term  competing  risk  problems.  Our  goal  is  to  employ  the  method  of 
competing  risks  to  model  the  fundamental  duel,  a  goal  we  will  attain  by  undertaking  three  tasks: 

1.  Express  the  time  to  system  failure  as  a  function  of  the  components. 

2.  Express  the  probability  that  one  component  fails  before  the  remaining  components. 

3.  Relate  these  expressions  to  the  fundamental  duel. 

We  will  address  the  first  two  tasks  as  we  introduce  the  general  methods  of  competing  risks,  and 
arrive  at  the  third  goal  as  we  close  the  section. 


Figure  5.  The  fundamental  duel  as  two  marksmen. 
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We  start  with  a  system  that  consists  of  k  >  1  components  without  imposing  a  relationship  among 
these  components.  We  further  assume  that  the  likelihood  of  two  or  more  component  failures  in 
the  same  instant  of  time  to  be  negligible.  David  and  Moeschberger  (1978)  state  that,  in  general, 
such  a  system  possesses  two  observable  properties:  (1)  the  time  when  failure  occurs,  and  (2) 
which  of  the  k  components  failed.  If  we  let  {L(.  }*=1  denote  the  unobservable  lifetimes  for  the 
respective  components,  we  can  then  define  the  system  lifetime  L  in  terms  of  these  lifetimes  by 
letting  L  =  min  ,  L2 , . . . ,  Lk } .  Undoubtedly,  if  L  exceeds  some  value  t,  then  each  L  must  also 

exceed  t;  thus,  David  and  Moeschberger  write 

P[L>t)  =  P(Ll  >t,L2>t,...,Lk  >t ) 

=  1-^(0  (31) 


where  F[  (t  )  is  often  termed  the  survivor  function.  David  and  Moeschberger  define  the 
conditional  failure  rate  for  the  system  as 


rL(t)=  f' (t)  =h^l. 
1  ~FL(t)  Fl(t) 


(32) 


They  further  interpret  rL(t)  to  mean  the  probability  that  a  system,  given  that  it  is  operational  at 
time  t,  will  fail  in  (t,t  +dt),  owing  to  a  failure  of  any  of  its  components.  One  can  also  define  the 
conditional  failure  rate  for  the  system  in  terms  of  the  failure  rate  for  each  component.  First,  let 
gi  (t)dt  for  i  =  1, 2,  •  •  • ,  k  denote  the  probability  that  a  system,  given  each  component  operational 
at  t,  fails  in  (t,t  +dt)  because  component  i  failed.  Accordingly,  David  and  Moeschberger 
determine  that 


=  (33) 

1=1 

Furthermore,  they  infer  from  eq  33  that  one  can  view  the  “health”  of  the  system  as  a  function  of 
the  health  of  each  component  (1978,  p.  4). 

We  are  now  ready  to  address  the  first  two  goals  we  set  at  the  beginning  of  this  section.  We  begin 
by  rewriting  eq  32  as 

rL(.t)  =  -^-lnF[(t),  (34) 

dt 


then,  from  David  and  Moeschberger,  we  determine  that 


Q(a,t)  =  exp 


(35) 
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Equation  35  defines  the  probability  that  the  system,  having  been  operational  at  a,  remains  so  at 
a  + 1.  With  this  result,  David  and  Moeschberger  subsequently  determine  the  probability  that  a 
system  fails  in  (a,t)  due  to  a  failure  of  component  i 

t 

Qi(a,t)  =  ^gi(T)Q(a,T)dr.  (36) 

a 

Equations  35  and  36  respectively  satisfy  the  requirements  of  our  first  two  tasks.  However,  we 
now  note  that  while  computing  rL(t )  for  a  general  system  is  often  clear-cut,  computing  gt(t)  for 
each  component  of  that  system  calls  for  some  knowledge  of  the  system. 

Until  this  point  in  our  discussion,  nothing  has  required  us  to  impose  a  relationship  upon  the 
component  lifetimes.  The  literature  on  competing  risks  encompasses  a  large  amount  of  work  on 
models  where  the  relationships  among  lifetimes  are  independent,  dependent,  or  some 
combination  of  both.  In  addition,  this  body  of  literature  also  provides  substantial  discussion  on 
the  problem  of  identifying  lifetime  models  with  particular  functional  relationships  among  the 
component  lifetimes.  We,  including  other  researchers  such  as  Williams  and  Ancker  (1963), 
assumed  the  relationship  between  a  marksman  and  his  target  to  be  a  passive  relationship.  In  this 
sense,  a  passive  target  is  one  that  is  unaffected  by  the  act  of  having  been  fired  upon,  but  as  Bathe 
and  McNaught  (1989)  observed,  this  assumption  is  questionable.  Still,  in  order  to  reach  our 
immediate  goal,  we  now  find  it  necessary  to  rely  upon  a  passive  target.  For  the  moment,  we 
shall  justify  this  assumption  by  noting  that  Tsiatis  (1975),  Peterson  (1976),  and  Miller  (1977) 
independently  suggest  that,  absent  specific  information  regarding  the  structure  of  the  system, 
without  an  assumption  of  independence  among  the  system  components,  we  cannot  determine  the 
required  joint  or  marginal  distributions  for  the  component  lifetimes.  However,  we  contend  that 
by  assuming  a  passive  target,  one  implies  a  relationship  between  a  marksman  and  his  target, 
which,  in  the  fundamental  duel,  naturally  leads  to  the  assumption  of  independence. 

Nevertheless,  as  we  conclude  this  paper,  we  shall  have  available  a  means  to  construct  a 
framework  that  will  allow  alternative  structural  relationships  between  the  combatants  in  a 
fundamental  duel;  thus,  we  believe  that  further  research  will  allow  us  to  relax  the  requirement  of 
independence. 

Therefore,  by  assuming  independence  among  the  subsystem  lifetimes,  eq  31  becomes 

P(L  >t)  =  Pm  >  t)P(L2  >*)•••  P(Lk  >  t ), 


or  equivalently 

FI:(t)  =  f\F;(t), 
1=1 

where  F?  (t )  =  1  -  Fi  (t)  =  P(Li  >  t).  We  also  have 


(37) 
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m 

F?(ty 


(38) 


gft)  = 


m 

Fl(t) 


uf;w=- 


j= i 
i*i 


where  fft)  is  the  density  of  Fft).  David  and  Moeschberger  (1978)  equate  the  last  quantity  in 
eq  38  to  r.(t) ,  and  term  it  the  cause  specific  failure  rate,  or  the  marginal  intensity  function. 

Thus,  for  a  system  with  independent  component  lifetimes,  they  establish  that 

g.(t)=n(t)  i  =  1,2,  -,k.  (39) 


David  and  Moeschberger  also  state  that  for  independent  component  lifetimes,  the  probability  of 
system  failure  because  a  given  component  fails  in  (t,t+dt)  is  the  same  whether  that  component 
is  one  of  k  or  the  only  component  present.  Consequently,  by  assuming  independence,  eqs  35  and 
36  respectively  become 


Q(a,t)  =  exp 


-}  2/4(r)<fr 

a  *=1 


=  exp 


*=1  a 


(40) 


and 

t 

Qi(a,t)  =  jri(T)Q(a,T)dT.  (41) 

a 

With  eq  41,  we  can  now  define  the  fundamental  duel  as  a  competing  risk.  Let  us  define  the 
events  R  and  B,  and  take  them  to  signify,  respectively,  that  Red  or  Blue  was  victorious  in  the 
combat.  Furthermore,  assume  the  probability  of  Red  and  Blue  simultaneously  rendering  each- 
other  mission  incapable  to  be  negligible;  thus  P(R )  +  P(B)  =  1.  We  signify  that  our 
fundamental  duel  is  a  one-on-one  engagement  when  we  set  k  =  2 ,  and  thus  from  eq  41,  we 
define 


P(B)  =  0,  (0,  oo)  =  Jr,  (T)  0(0,  r)dr 

:  (42) 

P(R)  =  02  ( 0,oo)  =  J  r2(r)Q(0,T)dr. 

o 

In  their  original  formulation  of  the  fundamental  duel,  Williams  and  Ancker  (1963)  set  the 
condition  for  victory  to  be  the  first  combatant  to  score  a  hit  upon  the  other.  We  can  model  this 
situation  by  computing  rft)  and  r2(t)  from  the  requisite  distributions.  The  simplest  method  is 
via  the  generalized  marksman  problem  (eq  24),  with  the  appropriate  data  for  Red  and  Blue; 
however,  we  could  also  set  pk  =  1  and  resort  to  the  damage  model  (eq  30)  as  an  alternative. 

Equation  42  completes  our  work  for  this  section;  however,  obtaining  useful  analytical  results 
with  a  competing  risk  formulation  can  prove  difficult.  In  the  following  section,  we  consider  two 
examples  that  Williams  and  Ancker  originally  studied  (1963).  We  shall  see  that  in  the  first 
example,  a  competing  risk  formulation  is  an  asset,  while  in  the  second  it  is  a  drawback. 

Although  in  the  latter  case,  we  present  an  alternative  approach  that  yields  a  solution  equivalent  to 
that  obtained  by  Williams  and  Ancker. 
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5.  Analytical  Validation  of  the  Renewal  Duel  Formulation 


In  the  previous  section,  we  constructed  the  fundamental  duel  as  a  competing  risk  problem;  we 
now  evaluate  that  approach  by  studying  the  same  examples  that  Williams  and  Ancker  considered 
(1963).  As  in  earlier  sections,  we  make  use  of  two  combatants,  Red  and  Blue,  engaged  in  a  one- 
on-one  combat;  we  declare  the  victor  to  be  the  first  combatant  to  score  a  hit  upon  the  other. 
Williams  and  Ancker  obtained  results  for  such  a  combat  that  correspond  to  the  probability  of 
Blue  surviving  engagement  with  Red  in  each  of  their  original  examples.  Our  goal  is  to  represent 
these  probabilities  by  employing  a  competing  risk  framework.  To  reach  this  goal,  we  require  the 
marginal  intensity  functions  for  Red  and  Blue.  In  order  to  find  these  functions,  however,  we 
must  know  how  the  interval  between  consecutive  hits  (not  shots)  is  distributed.  If  XB  and  XR 
describe  the  interval  between  consecutive  hits,  then  we  must  find  FB it)  =  P(XB  <  t)  or 
FR  (t)  =  P{XhR  <  t)  as  necessary.  However,  since  the  mathematical  models  describing  the  Red 
and  Blue  marksman  processes  are  the  same,  we  need  only  derive  the  marginal  intensity  function 
for  Blue,  and  then  to  obtain  Red’s  we  simply  replace  the  appropriate  subscripts.  In  the  remainder 
of  this  section,  we  employ  pR  and  pB  to  denote  the  respective  constant  probabilities  of  hit  for 
Red  and  Blue.  In  addition,  we  will  let  fB  (t)  and  fR  (t)  respectively  denote  the  densities  of 
FhB{t)mdFhR{t). 

We  first  consider  a  case  where  the  interval  between  consecutive  shots  taken  by  Blue  has  the 
exponential  probability  density  fB  it)  =  AB  exp(-/l Bt) ,  with  <pB  (s)  =  XB  /( s  +  A,B)  as  the  Laplace 
transform  of  fB(t).  Let  @B(s)  denote  the  Laplace  transform  of  fB(t).  We  obtain  @B(s)  by 
differentiating  eq  20  with  respect  to  t,  taking  the  Laplace  transform  of  the  result,  and  rearranging 
terms.  Thus,  we  write 


0»(s)=  Pb^X)  _  PbK 

1-(1  -Pb)<Pb(s)  s  +  pb  4' 


By  inverting  eq  43,  we  obtain 

f B  (0  =  Pub'll  ^Xp(— 


which  we  then  integrate  to  find 

t 

K  (t)  =  \  fB(t)dt  =  \-  exp  (~pBAB ). 


0 


(43) 


(44) 


(45) 
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We  can  now  obtain  the  marginal  intensity  function  for  Blue  by  turning  to  eq  38 


flit)  =  /?B/tBexp(-/ylB)  ^ 
l- Flit)  l-[l-exp i-pM]  B  B 


(46) 


We  find  rR  (t)  by  replacing  the  subscript  B  with  R  in  eq  46.  By  turning  to  eq  40,  we  compute  the 
time  until  combat  termination 


0(0,0  =  exp 


t  t 

-\rB{T)dr-\rR(T)dr 

0  0 


exp  [-(pA+ZV1*)']- 


(47) 


If  we  now  substitute  the  quantities  given  by  eqs  46  and  47  into  eq  42,  we  obtain  our  desired 
probability  of  survival 


P(B)  =  0,(0,  oo)  =  \rB(T)Q(0,T)dr 

0 


00 

=  \pbK  oxp[-(pBAB+pRAR)r]dT 

0 

_  Pb^b 
Pb^b  Pr^r 


(48) 


Thus,  we  find  that,  in  this  case,  a  competing  risk  framework  yields  the  same  result  as  Williams 
and  Ancker  (1963,  eq  6). 


Consider  the  second  example  of  Williams  and  Ancker,  a  case  where  the  interval  between 
consecutive  shots  taken  by  the  combatants  has  an  Erlang(2)  distribution.  We  describe  this 
distribution  and  its  associated  Laplace  transform  by 


fa  it)  =  exp(-ABt),  t>0,AB>0 


<PsBi*)  = 


s  4-  A, 


B  j 


(49) 


Recall  that  (P'lis)  denotes  the  Laplace  transform  of  fl it).  Once  again,  we  find  this  quantity  by 
differentiating  eq  20  with  respect  to  t,  taking  the  Laplace  transform  of  the  result,  and  rearranging 
terms.  Thus,  we  write 


Kis)  = 


pB<psBis) 


Pb^„ 


Pb^r 


(50) 


1-(1  -PB)<PsBis)  s2  +  2ABs  +  pBA2B  (s+agXs+bg)’ 

where  we  have  let  aB  =  AB  (1  +  -  pB ),  and  bB  =  AB  (1  -  -y/l  -  pB ).  By  inverting  eq  50,  we  obtain 


fhBit)  = 


'  pX  ' 

V  ~aB  J 


{exp(-«Bf )  -  exp(-bBt )} . 


(51) 
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We  then  integrate  eq  51  to  find 


t 

FhB(t)=\fhBm 


Pb*i 


2  A 


bB  uB  j 


exp  (~bBt)  exp  (~aBt)  1  1 


K 


a t 


=  1- 


B  B 

bB  exp {-aBt)-aB  exp (~bBt) 
bB -aB 


+  - 
bg  Clg 


We  now  obtain  the  marginal  intensity  function  for  Blue 

/*(') 


rB(t)  = 


1  -F*(0 

f  PjX  A 


{exp(-flBf)  -  exp(-6B#)} 


bB  exp(-aBf) -aB  exp {-bBt) 


bB  uB 


f  PjX  ^ 

\J)B  —UB  J 


{exp(-aBf )  -  exp(-bBt )} 


bB  exp (~aBt) -aB  exp (~bBt) 
bB-aB 

PgA2g  {exp(-aBt)-exp(-bBt)} 


(52) 


(53) 


bB  exp (~aBt) -aB  exp (~bBt) 

We  find  rR(t)  by  replacing  the  subscript  B  with  R  in  eq  53.  By  turning  to  eq  40,  we  now 
compute  the  time  until  combat  termination 


Q(0,t)  =  exp 


l  l 

-\rB(r)dT-\rR{T)dT 


(54) 


If  we  substituted  the  quantities  given  by  eqs  53  and  54  into  eq  42 

oo 

P{B)  =  \rB{r)Q{0,T)dr 


(55) 


and  then  carried  out  all  of  the  integrations,  we  would  obtain  the  probability  of  survival  we 
sought;  however,  that  integration  is  extremely  difficult.  In  this  case,  we  see  a  competing  risk 
formulation  that  leads  to  a  very  complicated  solution,  but  we  have  other  means  available  to  find 
an  analytical  solution  for  our  Erlang(2)  case. 
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Because  we  employed  the  criterion  that  the  first  combatant  to  score  a  hit  upon  the  other  becomes 
the  victor,  we  can  take  another  path  to  obtain  a  competing  risk  solution  to  the  Erlang(2)  case. 

Let  XhB  and  XhR  denote  the  respective  lifetimes  of  Blue  and  Red  respectively;  thus, 

Z  =  min  |  Xg ,  XhR  j  defines  the  length  of  the  combat.  If  we  let  the  indicator  function,  /,  take  on 
the  values 


J 1,  If  Blue  hits  Red  before  Red  hits  Blue, 
[O,  If  Red  hits  Blue  before  Blue  hits  Red, 


then  the  indicator  function  describes  the  victor.  Let  the  event,  B,  denote  that  Blue  is  victorious; 
thus,  we  can  define 

P(B)  =  P(I  =  l)  =  P(Xhg<XhR).  (56) 


We  have  available  an  expression  for  fR  it)  in  eq  51,  and  by  a  change  of  subscripts,  we  also  have 
available  an  expression  for  fR  (t)\ therefore,  we  can  immediately  compute  P(B)  by 


P(B)  =  P(Xhg<XhR  ) 

oo  oo 

=  \fhBi.t)\fhR{T)dTdt 

0  t 

1  1 


(dB  +  bR)bR  (bg  +bR)bR 


(tig  +  ttR  (bg  +  Ug 

1  1 
-  +  - 


(  pX  \ 

pXr  \ 

—dB  y 

V bR  ~aR  y 

(57) 


In  eq  57,  a  and  b  have  the  same  meanings  as  for  eq  51.  After  considerable  algebraic 
manipulation  eq  57  reduces  to  the  same  result  as  Williams  and  Ancker  obtained  for  the  Erlang(2) 
case  (1963,  eq  24). 

When  we  couple  this  section  with  the  previous  one,  we  conclude  that  we  can  construct  the 
fundamental  duel  using  a  competing  risk  framework.  However,  we  also  realize  that  a  competing 
risk  framework  does  not  necessarily  lead  to  the  nice  analytical  solutions  that  Williams  and 
Ancker  obtained  in  their  studies.  Consequently,  we  ask  ourselves  what  benefits  do  we  accrue  by 
employing  competing  risks  in  our  study  of  stochastic  duels,  and  we  respond  by  highlighting  the 
work  of  Ewing  et  al.  (2002)  who  modeled  the  complex  interactions  of  a  pest  with  a  California 
citrus  crop.  In  their  model,  Ewing  et  al.  assumed  a  hazard  function  that  was  reasonably  well 
behaved  over  time  (2002,  p.  41);  as  a  result,  they  described  the  occurrence  of  events  in  an 
ecosystem  via  this  hazard  function.  We  close  by  noting  that  while  a  competing  risk  framework 
may  not  offer  a  nice  analytical  solution,  the  framework  does  offer  the  potential  to  assist  in 
computer  simulation  studies. 
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6.  Stochastic  Duels  and  Layered  Survivability 


In  this  section,  we  introduce  the  concept  of  layered  survivability  (Horton,  1996)  in  order  to  build 
the  foundation  for  the  remainder  of  this  paper.  To  reach  this  goal,  we  undertake  two  tasks.  First, 
we  relate  both  the  generalized  marksman  problem  and  the  damage  process  to  the  encounter 
process  that  underpins  layered  survivability.  Second,  we  utilize  this  concept  as  a  tool  with  which 
to  modify  the  generalized  marksman  and  damage  processes  so  they  include  appropriate 
survivability  responses.  To  aid  our  discussion,  figure  6  reproduces  the  survivability  “onion” 
diagram  that  describes  layered  survivability.  We  note  that  layered  survivability  models  an 
encounter  between  combatants  who  can  interact,  but  whose  opportunity  for  interaction  is  not 
certain.  From  the  discussion  thus  far,  we  now  know  that  we  can  fill  the  role  of  encounter  with 
the  fundamental  duel  as  illustrated  by  figure  5. 

1.  Avoid  being  detected 

2.  Avoid  being  acquired 

3.  Avoid  being  hit 

4.  Avoid  being  perforated 

5.  Avoid  being  destroyed 


Figure  6.  The  layers  of  survivability  (Horton,  1996). 

First,  we  draw  attention  to  the  clear  and  the  filled  areas  of  figure  6.  Observe  that  the  boundary 
between  these  areas  represents  a  physical  break  in  an  encounter;  in  a  manner  of  speaking,  the 
boundary  represents  the  “skin”  of  the  platform.  If  one  thinks  of  the  vertical  dimension  in  figure 
6  as  a  measure  of  time  that  increases  as  we  move  closer  to  the  target,  then  above  the  fourth  layer, 
an  adversary  has  yet  to  hit;  hence,  we  may  only  consider  the  likelihood  of  that  event  by  modeling 
it  probabilistically.  However,  by  passing  into  the  fourth  layer,  the  adversary  has  hit  almost 
surely;  hence,  we  must  deal  with  the  damage  rendered  by  that  hit,  damage  that  can  range  from 
instant,  catastrophic  failure  to  no  damage  at  all.  Here,  we  model  events  before  the  fourth  layer 
via  the  generalized  marksman  process.  Necessarily,  the  outcome  of  the  generalized  marksman 
process  drives  our  need  to  evaluate  the  damage  a  target  sustains;  hence,  we  model  events  in  and 
beyond  the  fourth  ring  via  our  damage  model.  Accordingly,  figure  7  illustrates  how  these 
processes  are  coupled,  in  the  context  of  the  fundamental  duel  described  by  figure  5,  in  order  to 
capture  the  essence  of  layered  survivability;  the  processes  in  italics  represent  one  half  of  the 
encounter  and  the  non-italicized  processes  the  other  half. 
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Figure  7.  The  fundamental  duel  decomposed. 

For  purposes  of  this  section,  we  have  completed  our  first  task.  However,  as  we  discuss  figure  6, 
we  should  resist  the  temptation  to  characterize  each  layer  as  a  binary  event.  With  respect  to  a 
platform  and  its  crew,  only  two  layers  truly  represent  events  that  alter  the  crews’  perception  of 
the  platforms’  state  of  nature:  detection  and  destruction.  For  detection,  the  crew  either  knows 
with  certainty  that  there  is  an  adversary  in  proximity,  or  they  do  not.  If  the  crew  chooses  to 
engage  the  adversary,  they  either  know  with  certainty  that  they  have  defeated  the  adversary,  or 
they  have  not.  Thus,  in  the  fundamental  duel,  the  intervening  layers  are  the  exemplar  for  the 
wrestling  match  that  Clausewitz  described  in  his  quotation  regarding  war: 

War  is  nothing  but  a  duel  on  a  larger  scale.  Countless  duels  go  to  make  up  a  war, 
but  a  picture  of  it  as  a  whole  can  be  formed  by  imagining  a  pair  of  wrestlers. 

Each  tries  through  physical  force  to  compel  the  other  to  do  his  will;  his  immediate 
aim  is  to  throw  his  opponent  in  order  to  make  him  incapable  of  further  resistance. 

(1976,  p.  75) 

We  suggest  that  survivability  exhibits  a  dual  nature,  one  that  occurs  because  a  platform  both 
fights  and  defends  during  a  combat;  how  well  it  survived  was  a  function  of  successes  in  both 
facets.  Figure  7  illustrates  this  duality.  From  the  point  of  view  of  the  platform,  defense  responds 
to  the  arrow  directed  at  the  platform,  which  we  typically  refer  to  as  platform  survivability. 

Attack  is  represented  by  the  arrow  directed  at  the  opposing  platform,  which  is  platform  lethality. 

Now,  let  us  focus  on  platform  survivability  and  for  ease  of  discussion,  we  can  think  of  this  in 
terms  of  the  Red  generalized  marksman  process  coupled  to  a  Blue  damage  process.  The 
processes  of  detection  and  hit  are  intuitive  and  directly  considered  in  the  marksman  process; 
however,  we  must  choose  how  to  handle  acquisition.  In  some  cases,  we  may  model  acquisition 
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as  part  of  the  time  to  detect;  while  in  other  cases,  we  may  model  it  as  part  of  the  inter-firing  time. 
Consider  a  turreted  platform  that  has  detected  a  target.  In  this  case,  the  acquisition  time  would 
represent  the  time  to  slew  the  turret  to  aim  at  the  target;  hence,  we  would  model  the  time  to  slew 
the  turret  as  part  of  the  time  to  detect  and  fire  the  first  shot.  In  another  example,  consider  a 
platform  that  must  reacquire  after  each  shot.  In  this  case,  we  would  choose  to  model  the  time  to 
reacquire  as  part  of  the  inter-firing  time.  Now,  given  our  rudimentary  damage  model,  we  simply 
choose  to  consider  the  probability  of  lethal  damage  to  represent  the  effect  of  armor.  If  the  armor 
protects  us  from  damage,  we  survive  that  hit;  however,  if  the  armor  fails,  then  we  do  not  survive. 
A  more  detailed  damage  model  will  require  changes,  but  for  this  discussion,  this  notion  is 
sufficient.  Thus,  we  have  related  our  work  in  this  paper  to  the  concept  of  layered  survivability, 
and  completed  our  second  task.  Let  us  now  add  survivability  measures  to  our  renewal  duel 
formulation. 


7.  Stochastic  Duels:  Adding  Survivability  M  easures 


In  this  section,  we  again  study  a  one-on-one  engagement  between  a  combatant  and  his  target; 
however,  since  we  know  that  the  marksman  problem  of  Williams  and  Ancker  is  a  special  case  of 
our  generalized  marksman,  we  need  not  parallel  the  development  in  section  3.  First,  we  begin 
with  the  combatant’s  generalized  marksman  process,  and  add  survivability  measures  that  will 
disrupt  his  firing  process.  Next,  we  continue  with  the  target’s  damage  process,  and  add 
survivability  measures  that  perturb  the  interaction  of  the  combatant’s  projectile  and  the  targets 
armor.  Finally,  at  the  conclusion  of  this  section,  we  will  have  a  complete  renewal  duel 
formulation  with  survivability  with  which  to  model  the  fundamental  duel  of  Williams  and 
Ancker. 

7.1  The  Generalized  Marksman  with  Survivability 

We  begin  by  considering  a  combatant  as  he  fires  at  a  target.  We  have  modeled  this  firing 
process  via  the  generalized  marksman  problem  (eq  24).  Recall  that  the  combatant  begins  with  an 
unloaded  weapon  at  time  zero.  He  must  also  search  the  battlefield  to  find  a  target;  however,  now 
the  target  is  no  longer  passive,  as  we  will  allow  the  target  to  perturb  the  combatant’s  firing 
process.  We  let  {S-  ,I(Yi),I(Wi)}  denote  our  stochastic  process;  however,  we  augmented  the 
notation  of  section  3.2  by  introducing  W  to  reflect  that  the  process  also  depends  on  the 
survivability  response.  We  employ  the  indicator  function  in  the  following  sense 

I(W)  P’  the  survivability  measure  perturbed  the  shot, 

[0,  the  survivability  measure  failed  to  perturb  the  shot. 
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As  before,  we  employ  renewal  theory  to  model  this  process  in  which  shot  epochs  occur  {5/  }*=1, 
and  hits  occur  at  epochs  {Sk  Yk=r  X. ,  Xk,  and  F(  will  take  the  same  meanings  as  section  3.2. 
X.  is  distributed  as 


P(X*  <x)  = 


I  G'x(x), 


i  >  1, 
i=  1. 


Obviously,  the  shot  outcome,  F  ,  will  depend  on  the  survivability  response.  While  for  most 
practical  survivability  measures  one  should  not  assume  Yi  and  Xs.  to  be  stochastically 
independent,  we  will  start  by  making  this  assumption. 


We  begin  by  invoking  the  renewal  argument  that  M(t)  =  P(X’'  <  t ),  and  proceed  by 
conditioning  on  the  time  of  arrival  of  the  first  shot 


[Mm 


X  <t, 
X  >  t. 


(58) 


We,  again,  only  allow  for  two  shot  outcomes:  either  the  shot  hit  or  it  missed.  Thus,  we  describe 
the  state  of  the  process  after  the  first  shot  by 


{p(X*  <  f)l /(*!),  *'=*} 


l M(t-x),  /(Fj)  =  1, 

1  1,  /(F1)  =  0. 


(59) 


However,  now  process  outcome  also  depends  on  the  survivability  response. 

For  each  survivability  response,  we  allow  only  two  possible  outcomes:  success  or  failure.  Again, 
if  we  turn  to  Doob  (1994),  we  can  define  the  probability  of  a  successful  response  as 

Pw=  J  dFw,  (60) 

I(W,)= i 


where  Fw  is  the  distribution  function  for  the  employment  of  our  survivability  response.  For 
purposes  of  this  paper,  we  assume  that  pw  is  both  a  fixed  and  known  quantity.  Similarly,  we 
find  the  probabilities  of  hit  conditioned  on  the  success  or  failure  of  the  survivability  response  as 


P  Ills  J  dFy  |/(W.  )=1  , 
I(Yi)=l 

Puf  =  J  dFY 


Fl/(W;)=0- 


/(*;)= i 


(61) 


Again,  we  assume  that  both  pMs  and  ph[f  are  fixed  and  known  quantities.  We  can  combine  eqs 
60  and  61  to  find  the  probability  that  the  outcome  of  a  given  shot  was  a  hit  by  applying  the  law 
of  total  probability 

Ph*=PMSPw+Phif(l-Pj  (62) 
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where  phw  denotes  the  probability  of  hit  in  the  presence  of  survivability  measures.  Now,  eq58 
becomes,  after  applying  the  law  of  total  probability  over  the  range  shot  outcomes  and  combining 
the  result  with  eq  59 


1  lPh„+0-~Ph„)M(t-x), 

J"io, 


x  <t, 

X  >t. 


(63) 


Then,  by  applying  the  law  of  total  probability  over  the  distribution  of  shot  intervals,  we  obtain 

oo 

M(t)  =  J  \_Mit )  I  X[  =  x]dFsx  (x) 

0 

t 

=  J  [Phw  +  (1  -  Phw  )M(t  ~  x)]dF°  (x)  (64) 

0 

t 

=  Pufi'x  «)  +  (!-  P„.  x)dF‘x  (x). 

0 


We  will  leave  eq  64  as  the  principal  result  of  this  section;  however,  one  asks,  “what  about 
detection?”  We  address  that  question  via  an  example.  Wand  et  al.  (1993)  introduced  the  process 
of  detection  into  a  stochastic  tank  duel.  In  their  duel,  they  suggested  that  the  negative 
exponential  density  adequately  modeled  the  time  to  detect;  thus,  we  can  completely  characterize 
the  process  of  detection  by  one  parameter,  the  mean  time  to  detect.  If  we  assume  that  the 
distribution  Gsx  it)  in  eq  64  solely  represents  detection,  and  that  it  is  a  negative  exponential 
distribution,  then  the  survivability  of  our  platform  is  a  function  of  the  mean  time  to  detect. 

Hence,  increasing  survivability  implies  that  we  delay  detection  via  camouflage  or  other  means. 
Therefore,  without  any  further  development,  we  can  leave  eq  64  as  our  representation  of  the 
generalized  marksman  problem  with  survivability. 

7.2  A  Simple  Damage  M  odel  with  Survivability  M  easures 

We  now  consider  the  second  half  of  the  combatant-target  interaction,  a  process  we  modeled  via 
our  damage  process.  In  eq  30,  we  developed  a  damage  model  that,  in  a  rudimentary  fashion, 
incorporated  the  effects  of  armor.  In  this  model,  the  effect  of  the  armor  on  an  incoming 
projectile  was  all  or  nothing.  That  is,  if  the  projectile  perforated  the  armor,  then  the  target  was 
rendered  mission  incapable,  but,  if  the  projectile  did  not  perforate  the  armor,  then  the  target 
sustained  no  damage.  We  now  wish  to  model  survivability  measures  that  augment  armor 
protection  by  interfering  with  the  projectile,  for  example  an  active  protection  system.  Again, 
{Sk,I(Ak),I(Rk)}  denotes  our  stochastic  process;  however,  we  augment  the  notation  of 
section  3.3  by  introducing  R  to  reflect  that  the  process  also  depends  on  the  survivability 
response.  We  employ  the  indicator  function  in  the  following  sense 

f  1,  the  survivability  measure  interfered  with  the  shot, 

k  [O,  the  survivability  measure  did  not  interfere  with  the  shot. 
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For  this  process,  hit  epochs  occur  at  {5*  }*=1 ,  and  the  epochs  of  hits  that  render  a  target  mission 
incapable  occur  at  {Sd  }”=1.  Xk  and  Xd  will  take  the  same  meanings  as  section  3.3.  We  assume 
that  the  distribution  of  Xk  is  either  a  given,  or  is  obtained  from  the  generalized  marksman 
process  (eq  24  or  64  as  appropriate).  We  assume  that  I(Ak )  and  I (Rk )  are  independent  of  the 
intervals,  Xk,  and  set  the  goal  for  this  section  as  finding  Xf  in  terms  of  Xk  and  the 
survivability  response.  Because  we  assumed  that  the  survivability  response  and  the  effect  of  the 
armor,  to  reach  our  goal  we  need  only  find  an  expression  for  pk  in  eq  30. 


First,  we  define  some  probabilities.  Recall  from  section  3.3,  that  when  a  projectile  perforated  the 
target’s  armor,  it  caused  lethal  damage.  This  is  still  true;  however,  now  the  likelihood  that  a 
projectile  perforates  the  armor  is  altered  by  the  survivability  response.  In  section  3.3,  eq  30 
modeled  the  effect  of  armor,  and  equated  pk  to  armor  effectiveness.  We  now  alter  that  equation 
by  incorporating  the  survivability  response 


Pa  v  J  dl\\i  ( Rt  i-r 

/(A*)= 1 

P  \  /  =  J"  dF 'A\I(Rk)=0 , 

I(Ak)= 1 


(65) 


where  s  and /  denote,  respectively,  the  success  or  failure  of  the  survivability  response.  Given 
that  a  projectile  will  hit  almost  surely,  we  deploy  a  survivability  response;  a  response  that  is 
either  successful  or  not.  Thus,  we  define 

P„=  J  dFR.  (66) 

KRk)=i 


to  be  the  probability  that  the  survivability  response  successfully  interfered  with  the  incoming 
projectile.  In  eq  66,  FR  is  the  distribution  of  outcomes  for  the  survivability  response.  By 
combining  eqs  65  and  66,  we  can  now  compute  pk  as 

Pk  =PMsPR+PMf^-PR)-  (67) 


Now,  we  refer  to  section  3.3  and  immediately  write  the  solution  to  the  duel  with  armor  and 
damage 

t 

M(t)  =  pkFhx  (j f )  +  (1  -pk)\  M(t  -  x)  dFhx  {x)  (68) 

0 

where  we  invoked  the  renewal  argument  M(t )  =  P(X'‘  <  t).  Consequently,  eq  68  becomes  the 
primary  result  of  this  section,  and  we  note  that  Fx  it)  may  either  be  a  given  or  obtained  from  eq 
24,  or  its  analog  with  survivability  eq  64. 
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8.  Conclusions 


In  this  paper,  we  studied  a  fundamental  duel  composed  of  two  marksmen.  Initially,  we  turned  to 
renewal  theory  in  an  effort  to  reformulate  the  marksman  problem  of  Williams  and  Ancker 
(1963).  We  then  generalized  the  marksman  problem  to  include  the  process  of  target  detection 
developed  by  Wand  et  al.  (1993).  Finally,  we  developed  a  rudimentary  damage  model  that 
included  the  effect  of  platform  armor.  Thus,  having  established  the  pieces  of  a  renewal  duel 
model,  we  suggested  that  one  could  view  the  fundamental  duel  as  a  series  system  with  two 
components.  With  this  idea,  we  turned  to  David  and  Moeschberger  (1978)  and  created  the 
fundamental  duel  via  their  method  of  competing  risks.  Initially,  we  adopted  their  method  as  an 
alternative  means  of  computing  the  metric  we  desired:  the  probability  of  survival.  While  we 
achieved  the  immediate  goal  of  defining  the  probability  of  survival  in  terms  of  each  combatant, 
we  required  an  assumption  of  independence  between  combatants  that  we  hoped  to  avoid.  We 
justified  independence  with  the  work  of  Tsiatis  (1975),  Peterson  (1976)  and  Miller  (1977); 
however,  we  also  realized  that  need  for  this  assumption  was  implied  by  the  notion  of  a  passive 
target.  We  now  realize  that  the  assumption  of  a  passive  target  can  be  relaxed;  however,  the 
result  is  a  competing  risk  model  with  dependent  life  times.  While  later  work  will  address  this 
issue,  for  now  we  highlight  the  work  of  Langberg  et  al.  (1978).  Langberg  et  al.  studied  the 
problem  of  converting  a  model  with  dependent  lifetimes  into  an  independent  equivalent  while 
still  preserving  essential  qualities;  hence,  we  concluded  it  may  be  possible  to  overcome  the 
requirement  for  independence.  Finally,  we  returned  the  concept  of  layered  survivability  to 
establish  the  basis  for  the  remainder  of  the  paper.  There  we  related  our  generalized  marksman 
problem  and  the  damage  process  to  layered  survivability,  and  concluded  by  adding  survivability 
to  the  renewal  duel  model. 
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